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ABSTRACT 

A strong hard X-ray luminosity from a blazar flags the presence of a very powerful jet. If the 
jet power is in turn related to the mass accretion rate, the most luminous hard X-ray blazars 
should pinpoint the largest accretion rates, and thus the largest black hole masses. These ideas 
are confirmed by the Swift satellite observations of the blazar S5 0014+813, at the redshift 
z = 3.366. Swift detected this source with all its three instruments, from the optical to the 
hard X-rays. Through the construction of its spectral energy distribution we are confident 
that its optical-UV emission is thermal in origin. Associating it to the emission of a standard 
optically thick geometrically thin accretion disk, we find a black hole mass Af ^ 4 x IO^^A/q, 
radiating at 40% the Eddington value. The derived mass is among the largest ever found. 
Super-Eddington slim disks or thick disks with the presence of a collimating funnel can in 
principle reduce the black hole mass estimate, but tend to produce spectra bluer than observed. 

Key words: Quasars: general — radiation mechanisms: non-thermal — gamma-rays: theory 
— X-rays: general 



1 INTRODUCTION 

Most of the estimates of the black hole mass in quasars make use of 
the width of broad emission lines in their spectra and an empirical 
relation between the luminosity of the ionising continuum and the 
size of the broad line region (BLR; Kaspi et al. 2007, Bentz et al. 
2006, 2009). This, plus the assumption of virial velocities, allows 
to estimate the black hole mass (albeit with some uncertainties). In 
recent years the Sloan Digital Sky Survey (SDSS) of galaxies and 
quasars is providing the largest quasars samples suitable to study 
the black hole mass of quasars, and how their corresponding mass 
function evolves with redshift. We here use a different (and "more 
ancient") method to estimate the black hole mass, by directly inter- 
preting the optical-UV flux of a source as the emission produced by 
a standard accretion disk, namely a Shakura-Sunjaev (1973) multi- 
colour disk, emitting as a black-body at each annulus. 

We do this exercise for carefully selected sources: blazars that 
are very luminous in hard X-rays, above 10 keV. The rationale 
for this choice is the following. We know that the spectral en- 
ergy distribution (SED) of blazars is characterised by two broad 
humps, whose peak frequencies (in vF^) are a function of the ob- 
served bolometric luminosity of the blazars (the so called "blazar 
sequence"; Fossati et al. 1998; Ghisellini et al. 1998). Larger pow- 
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ers correspond to smaller peak frequencies and more dominance of 
the high energy peak to the low energy one. According to this sce- 
nario, the most powerful blazars should have their peaks in the far 
IR and in the ~1 MeV bands. The latter peak should also carry the 
bulk of the electromagnetic output. This has two important conse- 
quences: i) the hard X-ray luminosity, being close to the peak, is 
large; ii) the non-thermal (synchrotron) radiation of the first peak 
(being in the far IR) does not hide the accretion disk radiation peak- 
ing at optical-UV frequencies. Therefore in these kind of blazars 
we can study the radiation from the accretion disk directly (see e.g. 
Landt et al. 2008; Maraschi et al. 2008; Sambruna et al. 2007). The 
last and important steps link the observed luminosity to the power 
carried by the jet (in bulk motion of particles and fields) and the 
jet power to the mass accretion rate. The latter point, investigated 
by Rawlings & Saunders back in 1991, has been since then con- 
firmed by other groups and using different methods (e.g. Celotti et 
al. 1997; Cavaliere & D'Elia 2002; Maraschi & Tavecchio 2003; 
Sambruna et al. 2006; Allen et al. 2006; Celotti & Ghisellini 2008; 
Ghisellini & Tavecchio 2008, 2009; Kataoka et al. 2008). 

From the above it is clear that having a sample of distant FS- 
RQs detected in hard X-rays allows to pinpoint the most accret- 
ing systems, and thus large black hole masses and large accretion 
rates. The recently published list of blazars detected by the Burst 
Alert Telescope (BAT) onboard the Swift satellite (Ajello et al. 
2009, hereafter A09) has provided the first well constructed sample 
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suitable for our aim. It contains 26 Flat Spectrum Radio Quasars 
(FSRQs) and 12 BL Lacs detected during three years of survey in 
the 15-55 keV range. One of those, S5 0014+813 (z = 3.366), 
is exceptionally bright in the optical band. Adopting a cosmology 
with S7m = 0.3 and Q.a = /lo = 0.7 we have, in the optical, 
vLv ~ 10** erg s~^. This is the source discussed in this paper, 
with the aim to find the mass of its black hole and the correspond- 
ing accretion rate. 



2 THE BLAZAR S5 0014+813 

This FSRQ was discovered in the radio band by Kuhr et al. (1981), 
and it was soon noted as exceptionally luminous in the optical 
(Kuhr et al. 1983). The spectrum taken by Sargent et al. (1989) 
had a (extinction corrected) slope of a = 0.8 [Fiu) cx jy""] long- 
ward of the Lya line of equivalent width of 158 A. Fried (1992) 
find no excess of foreground galaxies in the direction of the blazar, 
so disfavouring the hypothesis of gravitational lensing. It is not po- 
larised in the optical (Kuhr et al. 1983), and showed very mild op- 
tical variability (max Am ~ 0.15 in 9 years; Kaspi et al. 2007). 
At a Galactic latitude is 18.8°, it suffers from a non negligible op- 
tical extinction. Schlegel et al. (1998) hsts E{B - V) = 0.19, 
corresponding to Av ~ 0.62. 

This blazars was discussed in detail by Bechtold et al. (1994), 
that also analysed ROSAT data and showed the overall SED. These 
authors estimated the black hole mass on the basis of the optical- 
UV luminosity, that was however severely underestimated. 

The results of the XMM-Newton observations on 23 Aug. 
2001 are presented by Page et al. (2005). An absorbed power-law 
of photon index F = 1.61 ± 0.02 best fitted the data, with a column 
jyhost ^ ^-^ g ^ Q -^g^ ^ ^q22 locatgjj at the redshift of the 

source, in addition to the Galactic one. The [0.3-10 keV] flux was 
Fx = 5.6 X lO^'^^ erg cm^^ s^^, corresponding to Lx ~ 6 x 10""^ 
erg s~^. Observed with VLBI, it showed no superluminal motion 
(Finer et al. 2007). 



3 SWIFT OBSERVATIONS AND ANALYSIS 

Swift observed S5 0014+813 in January 2007. We analyzed these 
data with the most recent software Swift_Rel3 . 2 released as 
part of the Heasof t v. 6 . 6 . 2 . The calibration database is that 
updated to April 10, 2009. 

The XRT data were processed with the standard procedures 
(XRTPIPELINE V . . 12 . 2). Source events were extracted in a 
circular region of aperture ~ 47", and background was estimated 
in a same sized circular region far from the source. Response matri- 
ces were created through the xrtmkar f task. We analysed the sin- 
gle observations separately and also summed them together. Each 
spectrum was analysed through XSPEC with an absorbed power- 
law with a fixed Galactic column density (iVf = 1.32 x 10^^ 
cm~^ from Kalberla et al. 2005) The computed errors represent 
the 90% confidence interval on the spectral parameters. The best fit 
photon index of the summed spectrum was F = 1.36 ± 0.11 with a 
~ 16 for 19 degrees of freedom and an observed (de-absorbed) 
flux Fx = 5.3 X 10"^^ (Fx = 6.1 x 10^^^) erg cm'^ s"\ Afl 
quantities are calculated in the [0.2-10 keV] band. 

We have also re-analysed the XMM-Newton observation us- 
ing the SAS software (v. 9 . . ). After screening the data 
for high background phases, we have a net exposure of 11 ks. 
Fitting the 0.2-10 keV energy range with a power law we obtain 



r = 1.44 ± 0.05, with N^""^ = 6 x 10^^ cm"^ (in addition to 
A'^^f'''), harder than obtained by Page et al. (2005). This is likely 
due to our better screening and the improved calibrations. A bro- 
ken power law (with A''h fixed to A^h ^') better fits the data (at the 
99.99% level), with a break at Ei, = 1.0 ± 0.2 keV and slopes 
Fi = 1.1 ± 0.2 and F2 = 1.43 ± 0.04 below and above E^, 
respectively (see Fig.[TJ. 

UVOT (Roming et al. 2005) source counts were extracted 
from a circular region 5"— sized centred on the source position, 
while the background was extracted from a larger circular nearby 
source-free region. Data were integrated with the uvotimsum 
task and then analysed by using the uvotsource task. 

The observed magnitudes have been dereddened according to 
the formulae by Cardelli et al. (1989) and converted into fluxes by 
using standard formulae and zero points from Poole et al. (2008). 
The source was detected in V, B and U (observed magnitudes 
V = 16.43 ± 0.05; B = 17.57 ± 0.06 and U = 18.2 ± 0.1, 
not corrected by extinction), while only upper limits were obtained 
in the remaining three filters (UVWl > 19.4; UVM2 > 19.6 
and UVW2 > 20.0, 3a hmits). 



4 THE SPECTRAL ENERGY DISTRIBUTION 

Fig.[T]shows the SED of S5 0014+813. The BAT data correspond 
to the average flux of the 3 years survey, while the UVOT and XRT 
data are the sum of three pointed observations. The grey empty 
symbols are archival data, while filled grey symbols and thick 
(optical) segments are from Bechtold et al. (1994). The magenta 
points are from IRAS (Moshir et al. 1991) and 2MASS (Cutri et 
al. 2003). The solid lines correspond to our modelling (see below). 
This source is not in the list of blazars detected in the first three 
months of Fermi (Abdo et al. 2009). We have estimated the cor- 
responding upper limit shown by the arrow. There must be a peak 
between the BAT and the Fermi energy range. There is a slight mis- 
match between the level of the BAT flux and the extrapolation from 
the XRT spectrum, while the XMM-Newton data well agree with 
the XRT ones. The BAT/XRT mismatch can be easily accounted for 
by considering that all blazars are very variable sources (despite the 
XKI /XMM-Newton coincidental resemblance) and the BAT flux is 
a 3-years average, while the shown XRT spectrum is the sum of 
three observations within three days. 

The optical-UV spectrum is dominated by a a narrow bump, 
that we interpret as the emission produced by the accretion disk. 
This is substantiated by three facts: i) the broad emission lines of 
the source are well visible (e.g. Sargent et al. 1989; Bechtold et al. 
1994; Osmer, Porter & Green 1994) and are therefore not swamped 
by the non-thermal continuum; ii) there is a general consensus to 
interpret the overall non-thermal SED of blazars, from the far IR 
to 7-rays, as due to a single population of electrons. If so, the rela- 
tively steep emission above the high energy peak is made by elec- 
trons that emit, by synchrotron, a correspondingly steep spectrum 
in the IR-optical-UV band. This leaves the accretion disk com- 
ponent unhidden, iii) A synchrotron and inverse Compton model 
that reproduces the observed fluxes in the UVOT optical-UV and 
~MeV bands fails to reproduce the very hard optical spectrum [that 
has a slope F{iy) oc (although the non-simultaneity of the IR- 
optical data leaves some uncertainties). 

Then we model the optical-UV flux by a standard multi- 
colour accretion disk, with a temperature profile given by (see e.g. 
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Figure 1. SED of S5 0014+813 together with the fitting models, with pa- 
rameters listed in Tab. [T] UVOT, XRT and BAT data are indicated by red 
symbols, while archival data (form NED) are in light grey. The magenta 
square symbols are IRAS and 2MASS data points. The dotted line is the 
emission from the IR torus, the accretion disk and its X-ray corona. The 
blue and orange lines are the sum of all components. The two models differ 
mainly for the location of the dissipation region of the jet: outside (blue) or 
inside (orange) the BLR. 



Figure 2. Zoom of the SED of S5 0014+813, The three dashed fines cor- 
respond to L(j = I/Edd for tfiree black hole masses (40, 20 and 10 billion 
solar masses, from top to bottom). The three solid lines correspond to = 
0.4LEdd for the same three black hole masses. The blue line correspond to 
Af = 4 X 10"' and = O.lTLEdd- Note that M = 4 X 10'^°Mq 
can reasonably account for the two states of S5 0014+813 observed in 1993 
(Bechtold et al. 1994; grey symbols and fine) and in Jan 2007 (red symbols). 
A still smaller accretion rate can (approximately) account for the 2MASS 
data. 



Frank King & Raine 2002): 



1/2- 



(1) 



where — rjMc? is the bolometric disk luminosity, Rs is 
the Schwarzschild radius and ctmb is the Maxwell-Boltzmann 
constant. The maximum temperature (and hence the peak of the 
disk vFv spectrum) occurs at i? ~ 5-Rs and scales as Tmax oc 
(Ld/I/Edd)^''''A/"^/*. The total optical-UV flux gives Ld [that of 
course scales as (Ld/iEdd) M]. Therefore we can derive both the 
black hole mass and the accretion rate. 

The results of changing both are shown in Fig. |2l where we 
show the accretion disk luminosity for a black hole of 10, 20 and 
40 billion solar masses accreting at the Eddington or at the 40% 
of the Eddington rate. With 7\f = 4 x 10^° Af© we can reproduce 
reasonably well both the current state observed by UVOT and the 
old data discussed by Bechtold et al. (1994|3. The two states would 
then differ because of the different accretion rate. 

In Fig. [T] we show two theoretical models, one accounting 
for the BAT data, but overproducing the ~10 keV flux, the other 
accounting for the entire XRT and XMM-Newton data range but 
under-producing the BAT flux. The two models, bracketing the 
XRT and the BAT states, are obtained with a minimal change of 
the values of the input parameters (listed in Tab.[TJ. The little dif- 
ferences have no impact on our conclusions. We used the model 
discussed in detail in Ghisellini & Tavecchio (2009), accounting for 
the presence of many sources of photons located externally to the 
jet. We would like to stress that the use of a particular jet model is 



^ Note that Bechtold et al. used a cosmology with Hq = 100 and go 
0.5, giving a smaller distance than used here. 



not crucial for our discussion. What is important is that the optical- 
UV emission is dominated by the accretion disk emission. Indeed, 
as Fig.[T]shows, the non-thermal emission from the jet contributes 
only for a few per cent of the total flux. 



5 DISCUSSION 

5.1 Comparisons with other black hole mass estimates 

The value of 4 x 10^° A/q for the black hole of S5 0014+813 is un- 
precedented for a radio-loud source. Among quasars studied with 
the velocity width and continuum luminosity method, there seems 
to a be a saturation value around 5 x 10^-10^^ Mq (Shen et al., 
2008), with very few black hole masses above W^'^Mq (see also 
Vestergaard & Osmer 2009; Vestergaard et al. 2008; Kelly, Vester- 
gaard & Fan 2008; Natarajan & Treister 2009). 

For one specific quasar, Q0I05-2634, Dietrich & Hamann 
(2004) estimated a black hole mass of (41.4 ± 12.2) x 10^° M@ 
using the line width and the luminosity/BLR size in Kaspi 
et al. (2000), but for the same object the estimate decreases to 
(5.2 ± 1) X 10^ Mq using the Mgll line and the McLure & Jarvis 
(2002) BLR size/luminosity relation. This testifies the large uncer- 
tainties related to this method of estimating the black hole masses. 

Our method is free from this kind of uncertainties, but relies 
mainly on three assumptions: i) the emission of the disk, apart from 
a cos 8 term, is isotropic; ii) a standard (optically thick, geomet- 
rically thin) disk and iii) a black-body emission at each radius. 
While the latter assumption is conservative (since the black-body 
is the best radiator, it is bound to give a lower limit to the derived 
masses and accretion rates), the first two hypotheses can seriously 
affect our mass determination. We will first discuss the implica- 
tions of having found such a large mass, assuming that 40 billions 



© 0000 RAS, MNRAS 000, 000-000 



4 G. Ghisellini et al. 





M 




P' 




B 


r 




7b 


Tmax 


Sl 


S2 


logPr 


logi^ 


logPe 


logPp 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


[13] 


[14] 


[15] 


[16] 


9.6e3 (800) 


4el0 


4.9e3 


0.1 


2.4e3 (0.4) 


0.22 


16 


3 


70 


3e3 


-1 


2.3 


46.3 


46.6 


45.7 


47.4 


4.8e3 (400) 


4el0 


4.9e3 


0.1 


2.4e3 (0.4) 


1.54 


11 


3 


60 


3e3 


-1 


3.6 


46.1 


47.4 


44.5 


47.1 



Table 1. Parameters used to model the SED. Col. [1]; dissipation distance in units of lO'^^ cm and (in parenthesis) in units of Rs ; Col. [2]: black hole mass in 
solar masses; Col. [3]: size of the BLR in units of 10^^ cm; Col. [4]: power injected in the blob (in the comoving frame), in units of 10*^ erg s"'^; Col. [5]: 
accretion disk luminosity in units of 10^^ erg s~^ and (in parenthesis) in units of iEddJ Col. [6]: magnetic field in Gauss; Col. [7]: bulk Lorentz factor at 
Rd\ss\ Col. [8]: viewing angle in degrees; Col. [9] and [10]: break and maximum random Lorentz factors of the injected electrons; Col. [1 1] and [12]: slopes 
of the injected electron distribution [(3(7)] below and above 7^; Col. [13]; [14]; [15] and [16]: jet power in the form of radiation, Poynting flux, bulk motion 
of electrons and protons (assuming one proton per emitting electron). These powers are derived quantities, not input parameters. For a detailed description of 
the parameters and the model see Ghisellini & Tavecchio (2009). 



solar masses is the real value, and then we will discuss how a non- 
standard disk can impact our mass estimate. 



5.2 Consequences of the huge black hole mass estimate 

The fact to have found such a huge black hole mass in a blazar has 
a simple and profound implication: since we selected it on the basis 
of the beamed non-thermal hard X-ray continuum, there must be 
many other sources like S5 0014+813 that are pointing in other 
directions, with a much fainter (de-beamed) non-thermal emission. 
The relative number scales as r^>100, where V is the bulk Lorentz 
factor of the jet. These misaligned sources are unnoticeable in the 
X-ray and 7-ray bands, but since the accretion disk is unbeamed 
(and i/its emission is isotropic), they should be detectable in optical 
all sky surveys, like the SDSS. Assuming that there are 100 sources 
as optically bright as S5 0014+813 in the 30,000 deg"^ of the sky 
excluding the Galactic plane with |6| < 15°, the surface density of 
these objects is of the order of E ~ 3.3 x 10"^(r/10)^ deg"^. We 
can estimate how many of these sources the SDSS can detect. For 
the part of the sky already monitored and covered by spectroscopy 
by the SDSS (i.e. 5,700 square degrees for quasars) we expect ~19 
sources like S5 0014+813. Having a few objects that are indeed 
that luminous, the SDSS results are in this respect borderline (e.g. 
Vestergaard 2009). 

We can also compare these estimates with the expectations of 
different models/correlations relating the black hole mass (MBH) 
with the velocity dispersion a, the mass of the host galaxy bulge 
and its dark mass halo. 

Assume first that the relations among these quantities are 
redshift-independent. Then, from Tremaine et al. (2002), a MBH 
of 40 billion Mq should correspond to cr = 824 km , and ac- 
cording to Ferrarese (2002) this yields Mhaio = 6.7 x lO^^M©. 
Adopting a Press & Schechter law we derive E ~ 0.07 deg~^. 

Alternatively, we can derive the mass of the host bulge by ap- 
plying the relation Mbh,8 ~ ^-^^Mh-a^cii proposed by Haring 
&Rix (2004), finding Mbuigo = 1.4 x IO^^Mq. This Mbuigo can 
be related to a by assuming the "fundamental plane of black holes" 
proposed by Hopldns et al. (2007), finding a ~ 10^ km s~^, and so 
Mhaio = 1.3 X 10^" Mq, corresponding to E ~ 7 x 10"" deg"^. 

The above estimates assumed that the correlations between 
black holes and their host are redshift-independent. However, the 
MBH—Mi,„igf. and the Mbh-<j correlations might evolve with the 
cosmic time, as suggested by McLure et al. (2006); Peng et al. 
(2004); Treu et al. (2004); Woo et al. (2006). IVIc Lure et al. (2006) 
suggested that the ratio Msj/ /Afbuigc can evolve as (1 + z)"^ . In 
this case, Mbuige ^ 2.3x 1O^^M0. Through the fundamental plane 
we find a ~1800 km s"^ and Mhaio = 5.1 x lO^M®. This large 



halo mass makes the surface density to drop to E ~ 2 x 10~ 
deg~^. Note however that such a high velocity dispersion is rather 
extreme, implying a very dense bulge, with a size of ~ 3 kpc, 
within a very large halo. 

Finally, Woo et al. (2008) proposed that, at a fixed velocity 
dispersion, the IVIBH scales with redshift as (1 -I- z)^. This means 
that, at large redshifts, we can find larger MBH within hosts with 
smaller bulge and halo masses. We then use the Tremaine et al. 
(2002) relation to find a = 257 km s~^, corresponding to Mhaio = 
3.5 X 10^^ Mq. This relatively small halo mass corresponds to a 
very large surface density: E 6000 deg"'^. 

We can conclude that the current predictions for the number 
density of the largest MBH are far from conclusive, essentially be- 
cause we are far in the exponential tail of the distributions, and 
small changes of the host properties can dramatically change the 
predicted numbers: while the first two estimates of E are not far 
from the limits derived from our (single) object, the last two are 
either too small or too large. 



5.3 Super-Eddington disks 

The observed large optical luminosity is at the base of our black 
hole mass estimate. While we can exclude that it is dominated by 
beamed jet emission, we discuss here the possibility that the under- 
lying accretion disk producing it might by non-standard. In the lit- 
erature, two main alternatives have already been proposed: geomet- 
rically thick and radiation supported disk, for super-Eddington ac- 
cretion rates (Jaroszynski, Abramowicz & Paczynski 1980) and the 
so-called slim disk, with accretion rates close to the Eddington one 
(Abramowicz et al. 1988). The latter can emit a super-Eddington 
luminosity because the advection of the flow helps gravity to sus- 
tain the radiation force, but are characterised by a relatively small 
height to radial distance ratio (H/R) that does not allow a strong 
collimation of the produced radiation. The foreseen emission is not 
a pure black-body, but a a modified one, since electron scattering 
is very important. This implies larger temperatures than in the stan- 
dard case (Szuszkiewicz et al. 1996), that might be a problem in 
our case, due to the relatively severe upper limits in the UV. 

In the thick disk case, instead, the emission is locally close to 
Eddington, but the presence of a narrow funnel can coUimate (via 
electron scattering) the radiation produced deep in the inner funnel 
into a narrow cone. Observing at small angles from the axis of the 
funnel, we can then have the impression of a super-Eddington lu- 
minosity (up to a factor of 10-20). However, as detailed in Madau 
(1988), it is the high frequency radiation that is boosted the most 
(being the one produced mostly in the inner funnel), and this can 
be a problem in our case. 
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6 CONCLUSIONS 

The optical-UV luminosity of the blazar S5 0014+813 exceeds 
IC** erg s^^ and through the construction of the overall SED it can 
be convincingly associated to the radiation produced by an accre- 
tion disk. The found black hole mass is 40 billion of solar masses 
accreting at the 40% of the Eddington rate, and in the past it might 
have reached the full Eddington rate. Since this source was found 
because of its relativistically beamed hard X-ray emission, there 
should be many other sources of same mass and accretion rate, but 
whose jet is pointing in other directions. Current theoretical esti- 
mates do not exclude this, but are very uncertain. 

There are ways to reduce the estimated black hole mass, in- 
voking non-standard accretion disks, that however tend to emit a 
spectrum bluer than observed. On the other hand, the real physical 
properties of these slim or thick disks may be somewhat different 
from the assumed ones, and we cannot rule them out. A more defi- 
nite mass estimate would also greatly benefit by simultaneous data, 
from the far IR to the UV. 

The fact that this very large black hole mass has been found in 
a radio-loud source may not be a coincidence, if the presence of a 
jet is a crucial ingredient for the transfer of the angular momentum 
of the accreting matter, allowing the black hole to grow faster, as 
pointed out by JoUey & Kuncic (2008). In this case the found black 
hole mass can be real, or else the presence of the jet induces a 
super-Eddington accretion rate, making the disk slim or helping the 
formation of a funnel. Either way, S5 0014+813 is an exceptional 
source, worth to be investigated further. 
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